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Abstract 

"Quantum trajectories" are solutions of stochastic differential equations. Such 
equations are called "Stochastic Schrodinger Equations" and describe random phe- 
nomena in continuous measurement theory of Open Quantum System. Many recent 
developments deal with the control theory in such model (optimization, monitoring, 
engineering...). In this article, stochastic models with control are mathematically 
and physically justified as limit of concrete discrete procedures called "Quantum 
Repeated Measurements". In particular, this gives a rigorous justification of the 
Poisson and diffusion approximation in quantum measurement theory with control. 
Furthermore we investigate some examples using control in quantum mechanics. 

Introduction 

One of the topic in Quantum Open System theory concerns the study of the evolution 
of a small quantum system Tio undergoing an indirect and continuous measurement (the 
small system is in contact with environment and the measurement is performed on the 
environment). In this context, the evolution of the system is usually described by classical 
stochastic differential equations called "Stochastic Schrodinger Equations". Essentially, 
two kind of equations are considered: 

1. The "diffusive equation" (Homodyne detection experiment) is given by 

dpt = L{pt)dt + [ptC + Cpt - Tr {pt{C + C)) pt]dWt (1) 
where Wt describes a one-dimensional Brownian motion. 
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2. The "jump equation" (Resonance fluorescence experiment) is 



dpt = L{pt)dt + 



JiPt) 



pt {dNt-Tr[J{pt)]dt) 



(2) 



Tr[J{pt)] 



where Nf 



J is a counting process with stochastic intensity Tr[J'{ps)]ds. 



Solutions of such equations are called "Quantum Trajectory"; they describe the evolu- 
tion of the state of the system perturbed by the continuous indirect measurement. 

Recent progress and developments in quantum optics and quantum information theory 
need a highest precision in experience process using measurement [IS] (sensitivity, minia- 
turization, optimization...). It imposed then quantum systems to be controlled. Two types 
of controls are usually considered in Quantum Mechanics: deterministic and stochastic 
control. 

A laser monitoring a qubit, i.e a two level atom system, is a basic example of determin- 
istic control application. The action of control is then characterized by the laser intensity 
(the term deterministic is related to the fact that we consider the intensity to be a deter- 
ministic function of time). Often, it is called "Open Loop Control". Such experiences are 
used in order to prepare systems in specific states for quantum computing. 

The notion of stochastic control is, here, directly connected with the procedure of 
measurement. Depending on the result of measurement, a control operation is performed 
in order to modify the evolution of the system. As a result of measurement is random in 
Quantum Mechanics (one of the axiom of the theory), the control becomes also random 
(it justifies the term of stochastic control). This kind of control is particularly used in 
engineering when some constraint of precision and optimization must be followed. Usually 
such control is called "Closed Loop Control" or "Feedback Control". 

From a theoretical point of view, an important question is to lay out a mathematical 
setup to modelize the action of control. The next step is to describe the evolution of 
controlled quantum system. 

Usually in the literature, in order to obtain and justify the classical stochastic Schrodinger 
equations ([T]) and ([2]), Quantum Filtering theory [U] or Instrumental Process theory [7] are 
used. Such techniques are based on the Hilbertian formalism of Quantum Mechanics and 
on the theory of Stochastic Quantum Calculus. It uses heavy analytic machinery and all 
the subtleties of the non commutative character of quantum probability (conditionnal ex- 
pectation in Von Neumann Algebra, partially observed system...). The starting point is the 
description of interaction between system and environment in terms of quantum stochastic 
differential equations (also called Hudson Parthasarathy Equations [23]). In order to apply 
such theory in the control setup, a theory in adequacy with the non commutative character 
have to be introduced. Even if it is satisfied, the derivation and the obtaining of stochastic 
Schrodinger equations with control is far from being obvious and intuitive (see [TT]) and 
there are less rigorous results. 

Recently, in the framework of the description of the interaction of a small system with 
environment (without measurement), in [1], the authors have introduced a discrete model 
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of interactions: "Quantum Repeated Interactions". The basic model is the one of a small 
system Tig in contact with an infinite chain of quantum system (S^jli Ones after the 
others each copy of Ti. interacts with Tio during a time h. 

Such approach of open quantum system yields a "good" and "useful" approximation 
model of continuous-time interaction models. Indeed by rescaling this interaction with 
respect to h, it is shown that models of interaction (described by stochastic quantum dif- 
ferential equations) can be obtained as continuous limits (h goes to zero) of discrete models. 
In the measurement setup, this approach has been adapted in [2S] and [2S]- In these arti- 
cles, it is then shown that classical quantum trajectories (solutions of equations and ([2])) 
can be obtained as continuous limits of discrete models of quantum measurement called 
"Quantum Repeated Measurements" (also called "Discrete Quantum Trajectory theory"). 
The idea of discrete indirect measurements consists in performing a measurement of an 
observable of H after each interaction between Hq and a copy of 7-^. A discrete quan- 
tum trajectory is then a discrete random process describing the evolution of the state of 
Tio undergoing such repeated measurements. In this case, the approach of the theory of 
stochastic Schrodinger equations via approximation results is essentially based on classical 
probability theory (there are no problem of commutativity). 

The main aim of this article is to adapt such technique in the framework of control. 
In this article, we present the notion of control in the model of quantum repeated inter- 
actions. In this setup, quantum repeated measurements give rise to a discrete models of 
indirect quantum measurement with control. By adapting convergence results of [25] and 
|2(j] . we obtain the description of stochastic Schrodinger equations with control. With 
this approach, all the problem concerning non commutativity are avoided and the physical 
justification of stochastic models is rigorous and intuitive. 

This article is structured as follows. 

The first section is devoted to present discrete models of quantum measurement with 
control theory. We remember the mathematical model of quantum repeated interactions. 
Next, we introduce an appropriate notion of control in this setup and by introducing 
the measurement principle, we obtain the description of discrete quantum trajectory with 
control. Next, in order to prepare final convergence results, we adapt and enlarge the 
asymptotic assumptions presented in [1] to the context of control. To investigate such 
problems, we focus on a central case in physical application: a two-level atom in contact 
with a spin chain. 

The second section is then devoted to continuous models. By applying the asymptotic 
assumptions on the two-level atom model of Section 1, we obtain two different discrete 
evolution equations (in asymptotic form) describing the evolution of the state of Hq. Each 
evolution equation describes the evolution of a discrete quantum trajectory with control 
for a specific observable. For each equation, we investigate the continuous limit equation 
and we show the convergence. 

In the last section, we present some applications of continuous models. In a first time, 
we investigate a model of a deterministic control: an atom monitored by a laser. By 
modehng a suitable interaction discrete model and by adapting the result of Section 2, 
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we obtain a continuous stochastic model for this concrete example. In a second time we 
present a use of stochastic control: "Optimal Control theory". We adapt classical results 
concerning this theory in the quantum language. 

1 Discrete Controlled Quantum Trajectories 

This section is devoted to the presentation of the model of discrete quantum trajectories in 
presence of external control. Here, we present a natural way for modeling control theory. 

1.1 Repeated Quantum Measurements with Control 

In a first time, let us remember the general context of quantum repeated interactions. 

A small system, represented by a Hilbert space Tio, is in contact with an infinite chain 
of identical independent quantum systems. Each piece of the environment is represented 
by Ti and interacts, one after the others, with Tio during a time interval of length h. For 
example, a copy of Ti. can represent an incoming photon or a measurement apparatus... 

The space describing the first interaction between Tig and Ti is defined by the tensor 
product TCq ® H. The evolution is given by a self-adjoint operator Htot on the tensor 
product. This operator is called the total Hamiltonian and its general form is 

Htot = Ho® I + I ® H + Hint 

where the operators Ho and H are the free Hamiltonians of each system. The operator 
Hint represents the Hamiltonian of interaction. This allows to define a unitary-operator 

U = e'^^'°\ 

and the evolution of states of Tio ® H, in the Schrodinger picture, is given by 

p ^ UpU\ 

After this first interaction, a second copy of 7i interacts with Tio in the same fashion and 
so on. 

As the chain is supposed to be infinite, the whole sequence of interactions is described 
by the state space: 

T = no®(^nk (3) 

fe>i 

where TCk denotes the fc-th copy of TC. The countable tensor product (S)fc>i '^fc means 
the following. Consider that 7i is of finite dimension and that {Xq, Xi, . . . , X„} is a fixed 
orthonormal basis of Ti. The orthogonal projector on CXq is denoted by |Xo)(Xo|. This 
is the ground state (or vacuum state) of H. The tensor product is taken with respect to 
Xq (for more details, see [1]). 
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The unitary evolution describing the k-th interaction is given hy Uk which acts as U 
on Ho <S> Ti.k, whereas it acts hke the identity operator on the other copies. If p is a state 
on r, the effect of the k-th interaction is then: 

p^UkpU^ 

Hence the result of the k first interactions is described by the operator 14 on ^{^) defined 
by the recursive formula: 

{Vk+i = Uk+iVk /,s 
Vo ^ I 

and the evolution of states is given by 

p^VkpV^. 

In this context, a main feature of this article is to present measurement and control 
theory. Let start by describing the control theory. An action of control consist in modifying 
the interaction at each new step depending on the previous step (this condition allows 
further to introduce stochastic control). Therefore if Uk is the unitary-operator describing 
the k-th interaction, it depends then on the time of interaction and on a parameter Uk-i 
which gives account of the control. Likewise this parameter depends on the interaction 
time; the operator Uk is then denoted by Uk{h,Uk-i{h)). 

The whole sequence u = {uk{h)) is called the "control strategy". In term of u, the k 
first interactions are then described by the unitary-operator V^: 

V,^ ^Ukih,Uk-iih))Uk-iih,Uk-2ih)) ...U^ihMh))- (5) 
Finally, the evolution in presence of control is given by 

p^v^ p vr. (6) 

Before to give a complete definition of control strategies, we have to introduce the repeated 
measurement model. 

Let us describe the basic procedure on each system of the chain. Let A be any observable 
on Tik with spectral decomposition A — X]j=o ^^ji consider its natural ampliation as an 
observable on F by: 

fe-i 

A*^ :=(g)/«)A® (g) / (7) 

j=0 j>k+l 

The accessible data are the eigenvalues of A'^ and the result of the observation is random. 
If p is any state on F, we observe Xj with probability 

P[to observe A^] = Tr[pP/], j G {0, . . . ,p}, 
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where the operator corresponds to the amphation of the eigenprojector Pj in the same 
way as ([7]). If we have observed the eigenvalue Aj, the "projection" postulate called "wave 
packet reduction" imposes the state after the measurement to be 

pkppk 

~ Tr[pPf]- 

Remark: This corresponds to the new reference state of our system. Another measure- 
ment of the same observable (with respect to this state) should give P[to observe \j\ = 
1. Hence only one measurement give a significant information; it justifies a principle of 
repeated interactions. 

Quantum repeated measurements with control are the combination of this previous 
principle and the successive interactions (EI). After each interaction, a quantum measure- 
ment induces a random modification of the state of the system. It defines then a discrete 
process which is called "discrete controlled quantum trajectory". The description is as 
follows. 

The initial state on F is chosen to be 

/i = P ® /^i 

where p is any state on TYq and each (3i = (3 where j3 is any state on Ti. The state after k 
interactions is denoted by p", we have: 

p" = /i vr. 

The probability space describing the experience is S^* where S = {0,...,p}. The 
integers i correspond to the indexes of the eigenvalues of A. We endow S^* with the 
cylinder cr-algebra generated by the cylinder sets: 

Ki,...,ik = ^ Vl^ /ui = ii, . . . ,uJk = ^fe}. 

Remarking that for all j, the unitary operator Uj commutes with all P'^ for all k < j. 
For any set {ii, . . . , ik}, we can define the following operator: 

/ifc (ii, ...,ik) = / ® ® . . . ® Pi, ® / . . . /i" / ® ® . . . ® Pj, (g) / . . . 

_ pk pi . u pi pk 

This is the non-normalized state corresponding to the successive observation of Aji, . . . , Aj^. 
The probability to observe these eigenvalues is 

P[to observe A^, • • • , A^fc] = Tr[/i^(ii, . . . , ik)]. 

By putting 

^'[Aii,...,iJ = P[to observe A^, . . . , Aj^], 
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it defines a probability measure on the cylinder sets of which satisfies the Kolmogorov 
Consistency Criterion. It defines then a unique probability measure on S^*. The discrete 
quantum trajectory with control strategy u on F is described by the following random 
sequence of states: 

p7: s^* B{r) 

From this description, the following result is obvious. 

Proposition 1 Let u be any strategy and (p") be the above random sequence of states we 
have for all u E T,^ : 



Tr 



Now at this stage, we can make precise the definition of control strategies which corre- 
spond to the case of deterministic or stochastic control mentioned in the Introduction. 

Definition 1 Let u = {uk{h)) be a control strategy and let (p") be a quantum trajectory. 

1. If there exists some function u from M to M" such that for all k : 

Uk{h) = u{kh), 

the control strategy is called deterministic. It is also called "open loop control". 

2. If there exists some function u from M x B{T) to M" such that for all k : 

Uk{h) = u{kh,p'^), 

the control strategy is called Markovian. It is also called "closed loop control" or 
'feedback control". If for all k we have Uk{h) = u{p^), this is an homogeneous 
Markovian strategy. 

The following theorem is an easy consequence of Proposition [1] and of the previous 
Definition. 

Theorem 1 For all control strategy u, the sequence (p")n is a non homogeneous Markov 
chain valued on the set of states ofT. It is described as follows: 

// p" = On then Pn+i takes one of the values: 

' Tr[Un+lih,UrmOnU:^,^,ih,Unih))Pl'+'y ' ' ' ' ' ' ^' 

With probability Tr [ f/„+i(/i, 0„ f/^+i(/i, ] . 

The discrete process (p") is called a controlled Markov chain. 
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Proof: Property of being a Markov chain comes from the fact that a control strategy 
is either deterministic or Markovian. For the two cases, the conclusion is obvious from the 
description of Proposition [TJ □ 

In general, one is only interested in the reduced state of the small system. This state 
is given by the partial trace operation. Let us recall what partial trace is. Let Z be any 
Hilbert space, the notation Trz [W^] corresponds to the trace of any trace-class operator W 
on Z. 

Definition-Theorem 1 Let Ti and /C he any Hilbert spaces. Let a be a state on the tensor 
product 7i ® /C. There exists a unique state rj onH which is characterized by the property: 

TrnlvX] = Trn^jc[a{X ® I)]. 

for all X G B{l-L). The state rj is called the partial trace of a onTi with respect to K,. 

For any state a on F, denote Eo[a] the partial trace of a on TYq with respect to 
®k>i ^fc- then define the discrete controlled quantum trajectory on TYq as follows. For 
allc^G T.^*: 

Pn(^) = MPli^)]- (8) 
Remark: We adapt Definition 1 by considering Markovian strategy defined on M x B{1-Lq). 
An immediate consequence of Theorem 1 is the following result. 

Theorem 2 For all control strategy u, the random sequence defined by formula ^ is a 
non-homogeneous controlled Markov chain with values in the set of states on 1-Lq. If = Xn 
then p"_^_^ takes one of the values: 



Er 



10 Pi Un+liK Unjh)) iXn ® P) U*+l{K Unjh)) I ® Pj 
Tr\ Un+l{h, Un{h)) {Xn ® /9) K+liK Un{h)) I ® P.,] 



0...p 



with probability Tr 



Un+i{h, Un{h)){xn ® (3)U*^^{h, Un{h)) Pi 



Remark: Let us stress that: 

{I®P^)U{Xn®l3) iI0P^) 

Tr[U{xn®P) U* (I^Pi)] 

is a state on Tio ® Ti.. In this situation, the notation Eq denotes the partial trace on TYq 
with respect to Ti. Moreover for each n, the operator f/„, which appears in the description 
of the transition of the Markov chain (p"), acts on T-Cq^H as the operators f/„ on Ho®Hn 

With the description of Theorem 2, we can express a discrete evolution equation de- 
scribing the discrete quantum trajectory (p"). By putting 



C'(p) 



En 



I^Pi Uk{h, Uk-i{h)) (p (g) p) jjl{h, Uk-i{h)) I® Pi 



0...P, 
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and li{uj) = li{u!k) for all uj G S^*, the discrete process (p") then satisfies 

for all G and all A; > 0. 

The following section is devoted to the deeply study of the equation ([H]) in a particular 
case of a two-level system in interaction with a spin chain. Next we come into the question 
of asymptotic assumptions. 

1.2 A Two-Level Atom 

The physical situation is described hy TCq = TC = C^. In this case, an observable A has 
two different eigenvalues: A = XqPq + XiPi (the case which only one eigenvalue is not 
interesting). The equation ([9]) is reduced to: 

P".iM = ^"''^"("^^ ^^M + ^"'"^(^"M) i^.i(^). (10) 

Pk+l %+l 

where = TrlC^'^'^^ {p^)] = 1 — Let now introduce the centered and normalized 
random variable 

We define the associated filtration on {0, 1}^: 

Tk = (T{Xi,i < k). 

So by construction we have E[Xfc+i/^fc] = and E[X|^^/JF^.] = 1. In terms of (X^) the 
discrete controlled quantum trajectory satisfies: 



+ 




Pk+l V 



Xk+1. (11) 



To give more sense to the equation ffTTj) . we have to express the terms £"''^'^^(p") in a 
more explicit way. For this, we introduce a specific basis. Let [Xq = Q, Xi = X) be an 
orthonormal basis of TCq = H = C^. For the space Hq ^H, we consider the following basis 

In this basis, the unitary operator can be written by blocks as a 2 x 2 matrix: 

Uk+i{h,Uk{h)) = 



Loo{kh,Ukih)) Lm{kh,Uk{h)) 
LiQ{kh,Uk{h)) Lii{kh,Uk{h)) 
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where each Lij{kh,Uk{h)) are operators on Hq. The reference state f3 of H is: 

(3 = \n){n\. 

The terms C^'^^^{p^) depends also on the expression of the eigenprojectors of the observ- 
able A. If the eigenprojector Pi is expressed as = [ ^-"^ ) in the basis (fi, X) of 

H, we have: 

^T^~^^iPk) = PmLooikh, Uk{h)) p" Ll^ikh, Uk{h)) + pl^Loo{kh, Uk{h)) p" L*Q{kh, Uk{h)) 

+ p\QLio{kh, Uk{h)) LlQ^kh, Uk{h)) + p\-^Lio{kh, Uk{h)) p^ LX^ikh, Uk{h)) 

(12) 

As the unitary evolution depends on the time length interaction h, the discrete quantum 
trajectory (p") depends on h. In Section 2, this dependence allows us to consider continuous 
time limit {h 0) of the discrete processes (p"). The next section is devoted to present 
the asymptotic ingredients necessary to obtain such convergence results. 

I. 3 Description of Asymptotic 

In this section, we present suitable asymptotic for the coefficients of the unitary operators 
Uk{h,Uk{h)) in order to have an effective continuous time limit from discrete quantum 
trajectories. Let h = 1/n be the length time of interaction, we have for (Uk) 

TT „ ( Loo{k/n,Uk{n)) Loi{k/n,Uk{n)) 

Uk^n^UkW) - y L,oik/n,Ukin)) L,,{k/n,Uk{n)) 

In our context, the choice of the coefficients Lij is an adaptation of the works of Attal- 
Pautrat in [1]. In their work, they consider only evolution of the type 



Uk+i{n) 



Loo{n) Lm{n) 
Lio(n) Lu(n) 



that is, homogeneous evolution without control. They have shown that 

V[nt] = U[nt]in) ...Ui{n) 

converges (in operator algebra) to a non-trivial process Vt (solution of a quantum stochastic 
differential equation), only if the coefficients Lij{n) obey certain normalization. In their 
case, these coefficients must be of the form 

i^oo(^) = ^ + ^ (-^^0 - ^^^') + ° (^) (13) 
Lio{n) = -^C + o(-], (14) 



n \n 
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where Hq is the Hamiltonian of Tio and C is any operator on C^. Hence, in the control 
context, the coefficients Lij{k/n,Uk{n)) must follow similar expressions. Let k be fixed, we 
put 



LoQ{k/n,Ukin)) = I + - ( -iHk{n,Uk{n)) - \Ckin,Uk{n))Cl{n,Uk{n))] +o(- 

n \ 2 J \n 

Loo{k/n,Uk{n)) = -^Ck{n,Uk{n)) + o ( - 



(15) 



where Hk{n,Uk{n)) is a self-adjoint operator and Ck{n,Uk{n)) is an operator on C^. It is 
straightforward that the expression ( fT3l) of Attal-Pautrat is a particular case of the previous 
expression. Finally, we suppose that there exist some function H and C such that 



H : 



i-^ X 



H2(C) and C : M+ x 
H{t,s) it,s) 



iVl2 

C{t,s 



where ]Hl2(C) designs the set of self-adjoint operators on and 

Hk{n,Uk{n)) = H{k/n,Uk{n)) 
Ck{n,Uk{n) = C{k/n,Uk{n)) 



(16) 



Furthermore we suppose that all the o are uniform in k. 

Now, we shall express the equation (fTTj) and (fT2|) with these asymptotic assumptions. 
As it was announced, we obtain two different behaviours depending of the choice of the 
observable. 



1. If the observable A is diagonal in the basis {Q,X), that is, it is of the form 
A = Xq { Q o)~'~'^Mo l)'^^ obtain the asymptotic for the probabilities 



-Tr 



n 



J{k/n,Uk{n)){pl{n)) 



+ o 



-Tr 



n 



J{k/n,Uk{n)){pl{n)) 



+ o 



n 



n 



The discrete equation (ITTl) becomes 

= ^L{k/n,Uk{n)){pl{n)) + o{-) 
n n 



+ 



J{k/n,Uk{n)){pl{n)) 



Tr 



J{k/n,Uk{n)){pl{n))) 



Pkin) + o(l) 



^/?"+i(^K+iW^fc+iH(17) 
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where for all state p, we have defined 

J{t,s){p) = C{t, s) p C*{t,s) a.nd 

L{t,s){p) = -z[H{t,s),p]-^{C{t,s)C*{t,s),p} + J{t,s){p). (18) 
2. If the observable A is non diagonal in the basis {Q, X), and if the eigenprojectors are 



express as Pq 



Poo Poi 
Pio Pn 



and Pi 



Qoo Qoi 
Qio qn 



we have 



Pk+i = Poo + p^ {poiC{k/n,Uk+i{n)) +piQC*{k/n,Uk{n))) 



+ -Tr 

n 



PkPoo {C{k/n,Ukin)) + C*{k/n,Ukin))) 



+ o I - 

n 



Qk+l 



qoo + —^Tr 

'n 



+ -Tr 

n 



Pk {qoiC{k/n,Uk{n)) + qioC''{k/n,Uk{n))) 



Pk Qoo {C{k/n,Uk{n)) + C*{k/n,Uk{n))) 
The discrete equation ( fTTl) becomes 



+ o I - I . 

n 



Pk+l Pk 



n 



n 



-L{k/n,Uk{n)){pl) + o ( - ) + e''C{k/n,uu{n))pl + e''' plC\k/n,Uk{n)) 



-Tr[pl {e''C{k/n,Uk{n)) + e-''C\k/n,Uk{n)))]pl + o{l) -^Xk+i{n) (19) 



where 6' is a real parameter. This parameter can be explicitly expressed with the 
coefficients of the eigenprojectors (Pj). By putting C6i(A;/?t,, Ufc(n)) = e^^C{k/n,Uk{n)) 
we have the same form for the equation ffTUl) for all 6, then we consider in the following 
that 6 = 0. The expression of L is the same as ffT^ . 

In order to prepare the final convergence result, in each case, we can define a process 
{P[nt]) which satisfies 

[ni]-l 

i=0 

[nt]-l 

= P0+ $^[C'+^(pr)+c^^'(pr)-pr] 

i=0 
[nt]-l 



i=0 



[nt]-l 



nt\~l 



j+1 



Po 



+ ]y{^/n,Ui{n),p'l)+ 2{t/n,Ui{n),p't)Xi+i 

i=0 i=0 



(20) 
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for some functions y and Z which depend on the description (fTTIl or (|T^ . 

Depending on the choice of observable, in the next section, we show that the process 
{P[nt]) converges to a solution of a particular stochastic differential equation. 



2 Convergence to Continuous Models 

In this section, starting from the description (1201) with a Markovian strategy and following 
the asymptotic (fTTl) and (fT9l) . we show that discrete processes (p[nt]) converge in distribution 
to solutions of stochastic differential equations. 

As in the classical case of stochastic differential equations, we show that the evolution 
of a quantum system undergoing a continuous measurement with control is either described 
by a diffusive evolution or by an evolution with jump. 

1. If (pt) denotes the state of a quantum system, the diffusive evolution is given by 



where (Wf) describes a one-dimensional Brownian motion. The function L is ex- 
pressed as ( ITS!) and 9 is defined by 



The functions L and are as (1181) . 

In a natural way, we call such equations "Controlled Stochastic Schrddinger Equations" 
and their solutions "Controlled Quantum Trajectories". 

For the moment, we do not speak about the regularity of the functions L, 9 and J'. 
This will be discussed when we deal with the question of existence and uniqueness of a 
solution for such equations. 

This question of existence and uniqueness is relatively important because this problem 
is not really treated in details in the literature. Moreover in the two cases, it is difficult to 
show that a solution of these equations is valued in the set of states (actually it is a essential 
property to solve the problem of existence and uniqueness). For physical considerations, 
this property is crucial otherwise these equations have no sense; we are going to see that this 
point can be deduced from the convergence result. Let us start by studying the diffusive 
case. 



dpt = L{t,uit, pt))ipt)dt + Q{t,u{t, pt)){pt)dWt 



(21) 



Q{t,a){p) = C{t,a)p + pC*{t,a) -Tr p\^C{t,a) + C*{t,a) j p 
for all t > 0, for all a in R and all operator p in M2(C). 



(22) 



2 
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2.1 Diffusive Equation with Control 



In this section, we justify the diffusive model 

dpt = L{t, u{t, pt)) {pt)dt + e (t, u{t, pt)) {pt)dWt 

of controlled stochastic Schrodinger equations by proving that the solution of equation (I2T1) 
is obtained from the limit of particular quantum trajectories {p[nt])- In the same time, we 
show that the equation (|2Tll admits a unique solution with values in the set of states. 

Start by investigating the problem of existence and uniqueness of a solution for fl2T|) . 
For the moment, let u be any measurable function which defines a Markovian strategy as 
it is expressed in Definition 2. Usual conditions concerning existence and uniqueness of a 
solution for SDE of type ( 1211) is that for all T > there exists a constant M{T) and K{T) 
such that the function L and G satisfy for all t < T and {p, p) G M2(C)^ : 

sup a)(/i) - a)(p)||, ||0(t, a)(^) - e(t, a)(p)||} <ir(T)||;t/-p|| 

sup {||L(t,a)(p)||, \m,a){p)\\} < M(T)(1 + ||p|| + ||a||). (24) 

Such conditions is called global Lipschitz conditions. However even in the homogeneous 
case without control, such conditions are not satisfied. Indeed, in the homogeneous situa- 
tion without control, for G we have 

e(t, a)(/i) = Q{p) =Cp + pC* - Tr [p{C + C*)] p. 

Such function is not Lipschitz. Nevertheless it is C°° and then local Lipschitz. Such 
property is used in the classical case to obtain the existence and the uniqueness of a 
solution for stochastic Schrodinger equations (see [25 and [2S])- In the non-homogeneous 
context with control, the local Lipschitz condition is expressed as follows. For all integer 
k > and all x G M, define the function (f)'' by 

The function (p^ is called a truncation function. Its extension on the set of operator on 
is given by 

4>\B) = (0^X^e(5,,)) + 20'=(/m(5,,)))o<.,,<i. 

Hence, the local Lipschitz condition for the functions L and 6 can be expressed as follows. 
For all T > and for all integer k > there exists a constant M^[T) and K^[T) such that 
the function L and G satisfy for all t < T and (/i,p) G M2(C)^ : 

\\L{t,a){^\p))-L{t,a){^\p))\\<K\T)\\p-p\\ 
\\Q{t,a){4>\p))-Q{t,a){4>\p))\\<K\T)\\p-p\\ 

sup {||L(t,a)(0'^(p))||, \m,a){^\p))\\} < M\T){1 + ||p|| + ||a||). (25) 
As a consequence, we have the following existence and uniqueness theorem. 
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Theorem 3 Let u he any measurable function. Let k > be an integer. Let {Q, J^, J-'t, P) 
be a probability space which supports a standard Brownian motion (Wt). Assume that L 
and satisfy the conditions (\25\\ . Let po be any 2x2 matrix. The stochastic differential 
equation 

Pt'' = Po+ f L{s,u{s,^'{p:'')){4>'{p:''))ds+ f e{s,u{s,4>'{p:''))){4>'{p:''))dW„ (26) 
Jo Jo 

admits a unique solution (p"'^). Furthermore the application t i— >• p"''^ is almost surely 
continuous. 

This theorem is just a consequence of the local Lipschitz condition ( l25i) (cf [29j). The 
process (p"''^) is called a truncated solution. The link between such solution and a solution 
of the equation (1^ without truncature is expressed as follows. Usually, we define the 
random stopping time 

Tfc = inf{t > 0/3(ij), i?e(p"'^(ij)) = k or Im{pt'''{ij)) = k} 

For any k > 1, we have > almost surely for po is a state and the almost surely 
continuity of (p"''^) (the coefficients of po satisfy namely \po{ij)\ < 1). Furthermore on 
[0, Tfc[ we have 

2k I" u,fc\ u,k 
<P [Pt ) = Pt ■ 

Therefore the process (p"' ) satisfies on [0,Tfc[ 

pT' = po+ t L{s,u{s,p:'')){p:''')ds+ te{s,u{s,p:'')){p:''')dw., (27) 

Jo Jo 

Hence the process (p"''^) solution of fpUl) is the unique solution of the equation fl?Il) on 
[0,Tfc[. 

In our situation, we will prove that = oo for all > 1 by proving that the process 
(p"''^) is valued in the set of states. Indeed if the process (p"'^) takes value in the set of 
states, we have for alH > 

2k f u,A;\ u,k 
<P {Pt ) = Pt , 

then Tk = oo. As a consequence the process (p"' ) satisfies for alH > the equation ( l2Ti) . 
The truncature method becomes actually not necessary, it just allow to exhibit a solution. 
As a consequence, we have to prove that the solution obtained with a truncature method 
takes value in the set of states. This property follows from the convergence theorem. 
Indeed, let assume that there is a discrete quantum trajectory (p"^]) which converges 

in distribution to (p"''^) (for some k > 1). Such convergence is denoted by 

u u,/c 

P[nt] =^ Pt ■ 

Therefore for all measurable functions V defined on M2(C), we have 

v(pr.,])^v(pr''^) 
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We apply it for the functions V(p) = Tr[p], for V(p) = p* — p and Vz{p) = {z,pz) for all 
2; G C^. By definition if p is a state we have from trace property Tr[p] = 1, from self- 
adjointness p* — p = and from positivity {z, pz) > for all z G C^. As discrete quantum 
trajectories take values in the set of states, these properties are then conserved at the limit. 
The limit process p"''^ takes then also values in the set of states. In |26], the problem of 
existence and uniqueness is proved independently of the convergence result. In this case, 
using convergence result is more practical because the equation is more complicated. Let 
us prove now the convergence result. 

Back to the description ( l20l) of discrete quantum trajectories, with asymptotic ( |T9l) in 
the case of a non-diagonal observable A and with a Markovian strategy, we have 



[ni]-l 



P[nt] 



Po 



k=l 



+ J2 - [L{k/nMk/n,p-,)){pl) + o{l) 



(28) 



+ 



[nt]-l 

E 

k=l 



[e {k/n, u{k/n, pH) (pD + -j^Xk^iin). 



From this description, we can define the following processes and functions: 



Un{t,W) 

e„(t, S) 

Ln{t, S) 



n 

k=l 

M 

n 

u{[nt]/n, W) 
Q{[nt]/n, s) 
L{[nt]/n, s) 



(29) 



for all t > 0, for all s G M and for all W G M2(C). 

By observing that these processes and these functions are piecewise constant, we can 
describe the discrete quantum trajectory (p"(t)) as a solution of the following stochastic 
differential equation 



Po + 



+ 



Po + 



+ 



Ln M„(S-, pl{S-)) {pl{S-)) + 0(1) dVn{s) 



e„ u„(.-, p:^ (.-)) {pi{s-)) + o(i) dw^{s) 



L„ ~^\p-{s-^)) {4>\pl{s-))) + 0(1 



e„ <p\pi{s-)) {<p\pi{s-))) + 0(1 



dVn{s) 



dWJs) 



(30) 
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for all k > 1. It appears essentially as a discrete version of equation fl2ip . This procedure 
was also used in [26j, but the equation without control is more simple and the convergence 
result needs less assumptions and arguments. 

Let us present the arguments in the control framework. In order to prove the conver- 
gence of this process to the solution of the equation (|26ll given by Theorem [3], we use a 
results of Kurtz and Protter ([20], [S]) concerning weak convergence of stochastic integrals. 
Let us fix some notations. 

For all T > we define ©[O, T] the space of cadlag process of M2(C) endowed with the 
Skorohod topology. 

Let Ti[0, oo) denote the set of non decreasing mapping A from [0, oo) to [0, oo) with 
A(0) = such that A(t + h) — X{t) < h for all t,h > 0. For any function G defined from 

M+ X M2(C) to M2(C), we define 

G : ©[0,00) X Ti[0,oo) — ^ ©[0,00) 
(X,A) ^ G(X)oA, 

such that for all t > we have G{X) o X{t) = G{X{t), Xx(t))- We consider the same 
definition for all other functions. We introduce the two following conditions concerning a 
function G and a sequence G„ as above. 

(CI) For each compact subset /C G ©[0, 00) x Ti[0, 00) and t > 0, 
sup sup ||G„(X,A)(s) -G'(X,A)(s)|| 0. 

(X,A) s<t 

{C2) For (X„, A„)„ G V[0, 00) x Ti[0, 00)/ sup ||X„(s) - X{s) \\ ^ 

s<T 

and sup |A„(s) — A(s)| for each t > implies 

s<t 

sup||G(X„,A„)(s)-G'(X,A)(s)|| ^0. (31) 

s<t 

Furthermore, recall that the square-bracket [X, X] is defined for a semi-martingale by 
the formula: ^ 

[X,X]t = Xf-2 [ Xs-dXs. 
Jo 

We shall denote by Tt{V) the total variation of a finite variation processes V on the interval 
[0,t\. The Theorem of Kurtz and Protter [21] that we use is the following. 

Theorem 4 Let {Hn,H) and {Kn,K) be two couple of functions which satisfy the condi- 
tions (CI) and (C2). Let (J^") be a filtration and let X^it) be a J^l^ -adapted process which 
satisfies 

Xnit) = X(0)+ / Hn{s,Xn{S-))dVn{s)+ [ K^S , X^S -))dWn{s) . (32) 

Jo Jo 
Let [Q, J-", J-'t, P) be a probability space. Let Xt be the unique solution of 

X„(t) = X(0)+ / H{s,Xs)ds+ [ K{s,X,)dWs (33) 
Jo Jo 
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where (Wt) is a standard Brownian motion on {Q, J-", J^t, P) . 

Suppose that (Wn, Vn) converges in distribution in the Skorohod topology to (W, V) where 
Vt = t for all t > and suppose 



sup < E 



[Wn,Wn] 



< oo, 



supjE" [Tt{Vn)]]<^. 



(34) 
(35) 



Hence the process (X„(t)) converges in distribution in V[0,T] for all T > to the 
process (Xt). 

We wish then to apply this theorem to obtain the convergence result for discrete quan- 
tum trajectories {p[nt]) described by ( l30l) . Concerning the convergence of the processes 
(Wn) and Vn we use the following theorem which is a generalization of Donsker Theorem 
(see 



Theorem 5 Let (M„) be a sequence of martingales. Suppose that 

lim E[\[M^,Mn]t-t\] =0. 

n— >oo 

Then Mn converges in distribution to a standard Brownian motion. 

In our context we have the following proposition. 
Proposition 2 Let {T^) be the filtration 

j-," = a(x„2< K]). 

The process {Wnif)) defined by is a J^^ -martingale. We have 

where iWt) is a standard Brownian motion. 
Moreover we have 



(36) 



supE 



[Wn^W, 



n\t 



< OO . 



Finally, we have the convergence in distribution for the process {Wn, Vn, ) to {W, V) 
when n goes to infinity. 

Proof: Thanks to the definition of the random variable Xi, we have E[Xj+i/J-!f] = 
which implies E 
property: 



i=[ns]+l ^H-'^t 

E[iy„(t)/^;] = w„(s) + E 



for t > s. Thus if t > s we have the martingale 

[nt] 



i=[ras]+l 



Wn{s\. 
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By definition of [Y, Y] for a stochastic process we have 

[Wn,W,,]t = W^itf-2 / Wn{S-)dWr,{s) = -y^Xf. 

-^0 ^ i=l 



Thus we have 



n ^ — ^ n ^-^ 

i=l i=l 

1 r 4-^ 



i=l 



Hence we have 



SW^^[[Wn,Wn]t] <t<00. 



Let us prove the convergence of (Wn) to (W). According to Theorem El we must prove 
that 

lim E[\[Mn,Mn]t~t\] =0. 

n— >oo 

Actually we prove a L2 convergence: 

lim E[|[M„,M„]t-t|2] =0, 

n— >oo 

which implies the Li convergence. In order to show this convergence, we use the following 
property 

B[Xf]=B[B[Xf/a{XiJ<t}]\=l 

and if i < j 

E[{Xf-l){X]-l)] = E[{Xf-l){X]-l)/a{X,J<j}]] 

= E[{Xf-l)]E[{X]-l) 
= 0. 



This gives 



E 



n t 



[nt] 



n 



[nt] 



^j:E[iXf-ir]+l^J2E[iXf-l)iX]-l 

i=l i<j 

[nt] 

i5:E[(.Y?-if], 



i=l 



Thanks to the fact that poo and goo are not equal to zero (because the observable A is not 
diagonal!) the terms E [{Xf — 1)^] are bounded uniformly in i so we have: 



lim E 

n—too 



[Wn,Wn]t- 



[nt] 



n 



0, 
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As ^ — > t in L2 we have the desired convergence. The convergence of (Wn, Vn) is then 
straightforward. □ 

In order to conclude to the convergence result by using Theorem HJof Kurtz and Protter, 
we have to verify conditions (CI) and (C2) for the functions appearing in the equation 
f lHUj) . We consider L„ defined by 

Ln{X) o (A)(t) = L„(A(t),n„(A(t),XA(t)))(XA(t)) + o(l) 

for all if: > 0, for all A G Ti[0, 00) and all cadlag process (Xt). Let us stress that in restriction 
to the processes {pt) which takes values in the set of states, the o are uniform in i^pt)-, we 
can then consider that the o are uniform for all processes. We define B„ in the same way. 
With this notation, we can express the convergence theorem. 

Theorem 6 Let J-"" be the filtration defined by (El). Let po be any state on Hq. Let (p"(t)) 
be the discrete quantum trajectory satisfying: 

Pnit) = P0+ f[L^{S-,Ur.{S-,p'^{S-m,p^{S-)) + o{l)]dV^{s) 

Jo 

+ [ [Qnis-,Unis-,p'^is-mp^is^)) + o{l)]dWn{s). (37) 
Jo 

Let k > 1 be any integer. Let (p"'^) be the unique solution of 

Pt'' = po+ fL{s,u{s,4>'{p:'')){4>\p:''))ds+ f e{s,u{s,4>'{p:''m4>'{p':''))dWs. (38) 

Jo Jo 

Assume the function u is sufficiently regular such that Ln, L and composed with 
the truncature function ip^ satisfy conditions (CI) and (C2). 

Then for all T > 0, the process (p"(t)) converges in distribution in V[0, T] to the process 
iPt)- 

Finally the process (p") is the unique solution of the controlled diffusive equation 

Pt=po+ f L{s,u{s,p:){p:)ds+ fe{s,u{s,p:'')){p:)dw,. (39) 

Jo Jo 

Proof: As the condition (CI) and (C2) are assumed to be satisfied, thanks to Propo- 
sition [2] and Theorem HI we have the convergence result. The final part of the theorem 
comes from the fact that the property of being a state is conserved by passage to the limit 
(see the remark at the beginning of this section). □ 

As regards conditions (CI) and (C2), the assumption for the function u is satisfied 
for example when u is continuous. By definition of the functions L„ and G„ conditions 
(Cl) and (C2) are namely satisfied for the functions L and satisfy the local Lipschitz 
conditions (1^ (used in Theorem [3] of existence and uniqueness) . 
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Hence, the model of diffusive stocliastic differential equation (12 ip for continuous mea- 
surement with control is physically justified by proving that solutions of such equations are 
obtained by limit of concrete discrete procedures. In the next section, we show a similar 
result by considering continuous limit of discrete quantum trajectories of type (ITTIl . 



2.2 Jump Equation with Control 



In this section, we investigate the convergence of a discrete quantum trajectory which 
comes from repeated measurements of a diagonal observable. 

In all this section, we fix a strategy u which defines a Markovian strategy. Furthermore, 
as in the diffusive case, we suppose that this strategy is continuous. Let A be any diagonal 
observable. With the use of description (|T7|) and (^U^, the discrete quantum trajectory 
satisfies 



[nt]-l 



P[nt] 



k=0 



+Tr[j{k/n,uik/n,p-,))ip-,)] + o(l 
J{k/n,uik/n,p-,)){p^) 



[nt]-l 
fc=0 



Tr[j{k/n,u{k/n,p^))ip^)] 



(40) 



Following the idea presented in article [25], we aim to show that the process {p[nt]] 
converges (n — >• oo) to a process (pt) which satisfies 



Pt 



Po + 



L(.-, u{s-, p:_)) {p:_) + Tr [J{s-, u{s-, p^)) (p^)] p^_ 



+ 



-j{s-,u{s-,p:_)){p:j ds 

J{s-,u{s-,p^_)){p^_) 
Tr[j{s-,u{s-,p-_)){p-_)] 



(41) 



p:- 



lo<x<Tr[J(s-,n(s-,p^_))(p^_)]^^(c?S, dx) 



where is a Poisson Point Process on M^. As a consequence, if the process (p") exists, it 



gives rise to the process (Nt) defined by 



'i-0<x<Tr[J{s-^,u(s-,p^_))ips-)]N{ds,dx) 



(42) 



which is a counting process with stochastic intensity 

fTr[J{s^,u{s-,p:J){p:_)]ds. 
Jo 
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In [23] , it is shown that the use of the Poisson process allows to give a mathematical 
sense to equation ([2SD and allows to define properly the process (Nt)- Actually in ([23]), 
the problem concerns the definition and the existence of the process {Nt). Indeed, the 
stochastic intensity of this process depends on the the solution of ( |23l) . In the way of 
writing the equation (1231) . in order to define A^^, it assumes implicitly the existence of the 
process solution whereas the equation is driven by (Nf). 

Now we consider the equation fBT]) as the jump-model of continuous time measurement 
with control. It will be justified later as limit of discrete quantum trajectories. 

For the moment, we deal with the problem existence and uniqueness of a solution for 
this equation. Let us denote 

R{t,a){p) = L{t,a){p)+Tr[J{t,a){p)]p-J{t,a){p) 

for alH > 0, for all a G M and all state p. It was obvious that (HTl) is equivalent to 

Pt = Po+ f R{s^,u{s^,p':_)){p:_)ds 
Jo 

+ 11 Q(s-,M(S-,p"_))(p"_)lo<x<Tr[J(s-,n(s-,p^_)){p._)]A^(c^S,rfx). 

Jo Jr 

In order to prove existence and uniqueness of a solution for such equations, a sufficient 
condition concerns the Lipschitz property for functions R an J'. In the same fashion of 
the diffusive case, this is not the case and a truncature method is used again. In the 
same way, we will have next to prove that the truncated solution takes values in the set of 
states. Regarding functions R an J', the conditions for the Poisson case are expressed in 
the following remark. 

Remark: As in the diffusive case, this remark concerns the regularity of the different 
functions. Firstly we suppose that R and J' satisfy the local Lipschitz condition fl25]) 
defined in Section 2.1. Secondly as the set of states is compact, we can suppose for the 
stochastic intensity that for all T > there exists a constant K{T) such that 

Tr[J{t,u{t,Xt)){Xt)]<K{T) 

for all t > T and for all cadlag process {Xt) with values in M2(C). This previous condi- 
tion implies the fact the stochastic intensity is bounded. Finally in order to consider the 
stochastic differential equation for all cadlag process, we consider the function 
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and the stochastic differential equation 

p-'' = po + f R{s-,u{s-,4>\p-L'))i4>'{p'::'))ds (44) 
Jo 

+ q{s-,u{s-,4>\p:l')){4>\p:l'))i,^.,^j,,^t,.^j^^^ 

where (p^ is a truncature function defined in Section 3.1. As in tlie diffusive case, if a 
solution of the equation (jHl) takes value in the set of states, it is a solution of the equation 
(14T|) . In addition to the diffusive case, we have to remark that if p is a state 

ReiTr[J{t,a){p)]) = Tr[J{t, a)ip)] > 

for alH > and for all a G M. 

Exactly in the same way as the diffusive case, if we show that a discrete quantum 
trajectory converges in distribution to a solution of the truncated equation (jHl), it involves 
that this solution takes values in the set of states. Let us first deal with the problem 
of existence and uniqueness of a solution for the equation (Hil) . We have the following 
theorem due to Jacod and Protter in 



Theorem 7 Let {Q,J-',P) be a probability space of a Poisson point Process N. The 
stochastic differential equation 



pT' = Po + / R{s-Ms-.~^\p:':^)){^\p':'-))ds (45) 
Jo 

+ Q{s-,u{S-,4>\p:L')){4>\p':L'))\^,^ne(TrlJis^Ms~,<P^ 





admits a unique solution p"''^ defined for alt>0. Furthermore the process 



ft 



O<x<_Re(rr[J(s-,«(s-,0fc(p^l'=)))(0*(p^i'=))]^(^'^' 

allows to define the filtration {J^t) where Tt = cr{Ns, s < t}. 
Hence the process 



Nt 



Re{Tr[J{s-, u{s^, 0^pr:')))(0'(p":')) 



ds 



is a Tt-martingale. 

The term (x)+ denotes max(0, x). Such theorem is treated in details in [25] for quantum 
trajectories without control. We give here a way to express the solution of (jS]) in a 
particular case. 
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Suppose that there exists a constant K such that: 



Re{J{t,u{t,Xt)){Xt)) 



<K, 



(46) 



for all if: > and all cadlag process (Xt). With this property we can consider only the 
points of contained in R x [0, K]. The random function 

M:t^iV(.,[0,t] X [0,K]) 

defines then a standard Poisson process with intensity K. Let T > 0, the Poisson Random 
Measure and the previous process generate on [0,T] a sequence {{Ti,^i),i G {!,... ,Aft)}}. 
Each Tj represents the jump time of the process (A/t). Moreover the random variables C,i 
are random uniform variables on [0,iir]. Let k > 1 he a fixed integer, we can write the 
solution of (jHl) in the following way: 



pT' = Po+ I Ris-,uis-,4>\p:L')m\p:L'))ds 



Aft 



+ Y,Q{n~Mn-A\Pr;-)){4>\Pr:')n 



0<i,<{Re{Tr[j{n-Mn-,4>HPr:-))MHp'^;'Lm+ 



i=l 



Aft 



<{Re{Tr[j(n-Mr^-,i^{p'',:t)mHfi'^:'L)]))+- 



(47) 



The general case is treated in details in [16]. Let us make more precise how the solution of 
dSj) is defined from the expression (H7I) in the particular case (H6l) . By applying Cauchy- 
Lipschitz Theorem, we consider the solution of the ordinary differential equation 



u,fc 



Pt 



P0+ / R{s-Ms-A\pT))W{p:'.))ds. 
Jo 



(48) 



It gives rise to the function 



Re{Tr[J{tMtA'{Pt''m4>'{Pt'')) 



Let define the first jump-time of the process (A^t). For this, we introduce the set 



Gt = {{x, y) e MVO < X <t,0 <y < 



Re{Tr[J{xM^,ct>\p-/)m\p-^')) 



and the random stopping time 



Ti = ini{t/N{Gt) = 1}. 



24 



As a consequence on [0,Ti[ the solution of (jHj) is given by the solution of the ordinary- 
differential equation fl48l) and p^^ is defined by 

U,fc U,k I r\lrr\ lrr\ U,fc \\l U,k \ 

Pt, = Pt^- + Q{Ti-,u{Ti-, p^^J){pj.^_}. 

We solve the ordinary differential equation after Ti with this initial condition, and by 
a similar reasoning we shall define a second jump-time. Thus, we construct a sequence 
of jump-time (T„). The boundness property (H6ll implies that the stochastic intensity is 
bounded. Hence, we can show limT„ = oo almost surely (see [16] or [25]). 

The solution of (jH]) is then given by the solution of the ordinary differential equation 

dp-'' = R{t,u{t,4>\pt'')m\Pt''))dt 

between the jump of the process A^^. The process A^^ corresponds to the number of point 
of the Poisson point process A^ included in the x axis and the curve 

ReiTr[Jit,uit,4>'{Pt'''m^\pt'')) 
The general case is more technical but can be expressed in the same way. 

In order to summarize the procedure to show the existence and uniqueness of solution, 
it is worth noticing that all the technical precaution are justified because we do not know 
that the equation preserves the property of being a state before showing the convergence 
result (it is the same problem in the diffusive case). 

Now we investigate the convergence result. In a first time, the way to proceed is the 
same as in ^25j. Next, we use another way to obtain the convergence result because we 
cannot applied the Theorem of Kurtz and Protter in this case. 

From expression (HOi) . define 

Pnit) = Ph 
[nt] 



Vnit) 



k=l 

[nt] 



n 

Rn{t,a){p) = R{[nt]/n,a){p), 
Qn{t,a){p) = Q{[nt]/n,a){p) 
Ur,it,W) = u{[nt]/n,W) 

for all t > 0, for all a G M and all W G M2(C). Hence the process (p"(t)) satisfies the 
stochastic differential equation 



Pnit) = / [R4s-,Ur,{s-,p:{S-)){p:{s-)) + o{l)]dVr,{s) 

Jo 

+ [ [Qn{s-,u^is-,p'^Xs^))ip^is^))+o{l)]dNr,{ 
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In this case, we do not have directly an equivalent of the Donsker Theorem for the 
process (A^„(t)). Because of the stochastic intensity of (Nt) which depends on the solution, 
it is actually not possible to prove the convergence of {Nn{t)) to {Nt) independently of the 
convergence of {p[nt]) to (pt). Hence we cannot applied the result of Kurtz and Protter 
used in the diffusive case. 

The convergence result is here obtained by using a random coupling method, that is, 
we realize the process {p[nt]) in the probability space of the Poisson Point Process A^. This 
method allows then to compare directly continuous and discrete quantum trajectories in 
the same probability space. It is described as follows. 

Remember that the random variables (if) satisfy: 

1^+1(0) = with probability p^+iin) = 1 - }^Tr[j{k/n,u{k/n, p^)){p^)] + o (1) 

1^+1(1) = 1 with probability qk+iin) = ^Tr[j{k/n,uik/n, p^)){p^)] + o (i) 

We define the following sequence of random variable which are defined on the set of 
states 

h+liV,^) = '^N{u,,Gkiri))>0 (49) 

where Gk{v) = < t < ^,0 < m < -n\n{Tr[C^o^\n){r])])}. Let po = P be any 

state and T > 0, we define the process (p^) for k < [riT] by the recursive formula 



+ 



(z>,+i(p-,.)-Tr[/:t+i(p-)]). (50) 



Thanks to properties of Poisson probability measure, the random variables (1^) and (z/^) 
have the same distribution. It involves the following property. 

Proposition 3 Let T be fixed. The discrete process {Pk)k<[nT] defined by {\5(J^ have the 
same distribution of the discrete quantum trajectory {p^)k<[nT] defined by the quantum 
repeated measurement. 

The convergence result is then expressed as follows. 

Theorem 8 Let T > 0. Let {Q, JF, P) be a probability space of a Poisson Point process N. 
Let {p^nt])o<t<T discrete quantum trajectory defined by the recursive formula i\5U^ 

Hence, for all T > the process {p^nt])o<t<T converges in distribution in V[0,T] (for 
the Skorohod topology) to the process (pf) solution of the stochastic differential equation: 

Pt = P0+ [ R{s-,u{s-,pt){p:_)ds (51) 

+ / / Qis~,u{S-, p'^_){p'^_)lo<x<Tr[J(s-,u{s-,P^-)){p^-)])Ni-^ ds, dx). 



26 



This theorem rehes on the fact that the process {p^nt]) satisfies the same asymptotic of 
the discrete quantum trajectory (pj^^]); we have namely 

[nt]-l ^ 

Pint] = Po+ J2 -[R{k/n,u{k/n,p'^)){p'^) + o{l)] 

k=0 
[nt]-l 

+ J2 [Q{k/n,u{k/n,pm~Pk) + °{^)]h+i{Pk,-)- (52) 

k=0 

The complete proof of this theorem is very technical. The idea is to compare the 
discrete process {p[nt]) with an Euler Scheme of the solution of the jump-equation. More 
details for such techniques can be found in [25] where the case without control is entirely 
developed. 

In the next section, we expose examples and applications of such stochastic models. 

3 Examples and Applications 

This section is devoted to some applications of quantum measurement with control. In a 
first time, by a discrete model, we justify a stochastic model for the experience of Resonance 
fiuorescence. The setup is the one of a laser driving an atom in presence of a photon 
counter. In a second time, we present general results in Stochastic Control theory applied 
to quantum trajectories. 

3.1 Laser Monitoring Atom: Resonance Fluorescence 

We here describe a discrete model of an atom monitored by a laser. A measurement is 
performed by a photon counter which detects the photon emission. The setup of repeated 
quantum interactions is described as follows. 

The length time of interaction is chosen tohe h = 1/n. Let us describe one interaction. 
Here we need three basis spaces. The atom system is represented by Hq equipped with a 
state p. The laser is representing by {li}, p}) and the photon counter by (H'^,P'^). Each 
Hilbert space are endowed with the orthonormal basis {Q, X) and the unitary operator 
is denoted by U. The compound system after interaction is: 

and the state after interaction is: 

a = U{p® p^ ® (3')U\ 

The appropriate orthonormal basis Tio^H} ^l-f^ , in this case, isfi(8>fi(8>f2, X®VL®VL, VL® 
X (g) fi, X (g) X (g) fi, (g) (g) X, X (g) (g) X, n ® X (g) X, X (g) X (g) X. As in the presentation 
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of the discrete two level atom in contact with a spin chain, the unitary operator is here 
considered as a 4 x 4 matrix U = (-^^i,j('^))o<i,i<3 where each Lij{n) are operator on TIq. 
If the different state of the laser and the counter as of the form 

for the state a = {auv)o<u,v<3 after interaction, we have 

ttuv = (aLuo{n)p + hLui{n)p^LlQ{n) + {cLuo{n)p + dLui{n)p^Ll^{n). (53) 

The measurement is performed on the counter photon side. Let A denotes any observ- 
able of Ti^ then I ® I ® A denotes the corresponding observable on 71^ ® H} ® l-C^ . We 
perform a measurement and by partial trace operation with respect to li} ® Ti'^ we obtain 
a new state on Tig- 

The control is rendered by the modification at each interaction of the intensity of the 
laser. This modification is here taken into account by the reference state of the laser. 
The reference state at the fc-th interaction is denoted by In the continuous case of 
Resonance fluorescence, the state of a laser is usually described by a coherent vector on 
a Fock space (see [9]). From works of Attal and Pautrat in approximation of Fock space 
(PP,[21]), in our context the discrete state of the laser can be described by 

, ^ / a{k/n) b{k/n) \ ^ 1 ( _ 1 h{k/n) \ 

\c{k/n) d{k/n) J l + \h{k/n)\^\h{k/n) \h{k/n)\^ J ' ^""^^ 

The function h represents the evolution of the intensity of the laser and depends naturally 
on n. 

If pk denotes the state on Tio after k first measurement, the state 

a>'+\n) = («L+'(^))o<,,,<3 = f/fc+i(ri)(p. ® pi ® P'WUM 
after interaction satisfies 

= (a{k/n)Luo{n)pk + b{k/n)Lui{n)pk^Llf^{n) 
c{k/n)Luo{n)pk + d{k/n)Lui{n)pk)Ll^{n) 



Remark: Let us stress that is not directly the framework of Section 1. Here, the 
control is namely not rendered by the modification of the unitary evolution. Moreover 
the interacting system is described by iJH} ® 'H'^,pi ® P) and p^ <^ P is not of the form 
|Xo)(Xo| as in Section 1. In order to translate this setting in the case of discrete models 
of Section 1, one can use the G.N.S Representation theory of a finite dimensional Hilbert 
space ([I9j,[l8]). This theory allows to consider the state /z^ ® /? as a state of the form 
|Xo)(Xo| in a biggest Hilbert space. The G.N.S representation modifies then the expression 
of operator Uk, and the control expressed in is again expressed in the new expression 
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of Uk (see [2] for more details). In our case, we do not use such theory because it is more 
exphcit to make directly computations to reach the discrete equation in asymptotic form. 



Let us present the result. The principle of measurement is the same as in Section 1. 
The counting case is also given by a diagonal observable of Tic- We shall focus on this 
case which renders the emission of photon ([9]). The asymptotic for the unitary operator 
follows the asymptotic of Attal-Pautrat in [3] . Let 6ij = 1 if z = j we denote: 

The coefficients must follow the convergence condition: 

lim n^''{Lij{n) - 6ijl) = Uj 

n— »oo 

where Lij are operator on Tio- 

Let Pq = ^"^^ Pi = be eigenprojectors of a diagonal observable A. If 

denotes the random state after k measurements we denote: 

= alt\n) + alt\n) (55) 

This is namely the two non normalized state, the operator £Q^^(pfc) appears with proba- 
bility pfc+i = Tr[4+i(pfc)] and C\-^\pk) with probability q^+i = Tr[4+i(pfe)]. 

From works of Attal-Pautrat in approximation and asymptotic in Fock space, we put 

h{k/n) = -^fik/n) + o (- 



where / is a function from M to C. In the same way of Section 2, we assume that the 
intensity of the laser / is continuous. 

With the same arguments of Section 1, the evolution of the discrete quantum trajectory 
is described by 



Pk= + 

Pk+i 



1^+^ (56) 



Pk+l Qk+l 

For a further use, convergence results will be established in the case Loi = —L\q, and 
Lii = L21 = L31 = -L30 = 0. Conditions about asymptotic of U and the fact that it is a 
unitary-operator imply that 



Lo, = -{iH + \Y.L^oUo) (57) 

i=l 

In the same way of Section 2.2, convergence result in this situation is expressed as follows. 
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Proposition 4 Let {Q, J-", J^t, P) be a probability space of a Poisson point process N on 



The discrete quantum trajectory {p[nt])o<t<T defined by the equation (E^ weakly con- 
verges in 'P([0, T]) for all T to the solution of the following stochastic differential equation: 



Pt 



Po 



/■* r 1 ^ 



10 



+ [f{s-)LioPs^ - f{s-)L*Q, p,_] - Tr[L2oPs-L2o]Ps 



ds 



+ 



-ps- + 



L20PS-L 



20 



Lo<x<Tr[L2oPs-L^ 



rr[L20Ps-^20JJ 

Proof: For example we have the following asymptotic for £Q^^(pfc): 
^o{Pk) = 



N{dx, ds). 



(5J 



1 

Pk + - 



n 



k — k 

Loop + pL'oo + LiopLlo + f{-)[LoiP + pLU + /(-)[^ioP + P^oil 



o I - 

n 



This above asymptotic, the condition about the operator Ljj and the theorem ([H]) prove 
the proposition. □ 

The stochastic differential equation ( l58il is then the continuous time stochastic model of 
Resonance fluorescence. In this model, the control is deterministic. Before to give an appli- 
cation of stochastic control , let us briefly expose a use of the laser monitoring atom model. 

Consider the special case, where the Hamiltonian H = 0. Let put 



C 



1 




10 — kiC, L20 — Kc 



kfC, 



with + l/ccP = 1- The constant kf and kc are called decay rates ([9]). 

Without control, the stochastic model of a two level atom in presence of a photon 
counter ([25]) is given by: 



Po + 



+ 



L 



1 



{C, Ps-} + Cps-C' - Tr[Cps-C*]ps 



ds 



-Ps- + 



Tr[Cps-C* 



'^Q<x<Tr[Ci,s-C*]N{dx, ds). 



(59) 



Let denote A^^^ = j'l ^o<x<Tr[CtJ.a-C*]N{dx, ds) and T = inf{t > 0; A^t > 0}. In ^ it 
was proved that: 

'10 




Pt 



\n){n\. 



(60) 
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for all t > T. Physically, it means that at most one photon appears on the photon counter. 
The mathematical reason is that if we write the equation fl59p in the following way: 

Jo Jo 

we have for fi = \Q){Q\ 

= vI/(^) = 0. 

The state is an equilibrium state. 

In the presence of laser, the control / gives rise to the term [fLiQ — fL^Q, .] = [hfC — 
kifC*, .]. Hence if = \Q){Q\,we still have = but we do not have anymore \E'(/x) = 
and the property (1601) is not satisfied. The state is no more an equilibrium state. 

As a consequence it is possible to observe more than one photon in the photon counter. 

In the next section we deal with general strategy and the particular problem of optimal 
control. Considerations about optimal control is an interesting mean to point out the 
importance of Markovian strategy. 



3.2 Optimal Control 

This section is then devoted to what is called the "optimal control" problem. It deals with 
finding a particular control strategy which must satisfy optimization constraint. In this 
section, we give the classical mathematical description of such problem and investigate gen- 
eral results in the discrete and in the continuous model of controlled quantum trajectories. 
Let us begin with the discrete model. 



3.2.1 The Discrete Case 

We come back to the description of a discrete quantum trajectory for a two-level system 
as a Markov chain. 

Let n be fixed, thanks to Theorem 2, a discrete controlled quantum trajectory (p") is 
described as follows. Let p be any state, if p" = p then p"^^ takes one of the values: 

^u,fc/ ^ ^ L,o{k/n,Uk{n)){p)L*Q{k/n,Uk{n)) -^^^ 
' Tr[Lio{k/n,Uk{n)){p)L*Q{k/n,Uk{n))] 

with probability, 

Pk+iip) = Tr[LoQ{k / n, Uk{n)) (p) L^Q^k / n, Uk{n))] foi i = 
Qk+iip) = Tr[LiQ{k/n,Ukin)){p)LlQ{k/n,Uk{n))] hi i = 1. 

With this previous description, the property of a strategy (uk) can be enlarged. We 
can namely consider more general strategies such that for all k the term Uk depend on all 
(pi) for i < k. We define U the set of all admissible strategies which satisfy this condition. 
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Let us stress that in this situation, the discrete quantum trajectory is no more a Markov 
chain because the strategy at time k depends on all the past of the strategy. 

With this remark concerning the definition of strategies, we can expose the general 
problem of "optimal control". In this article, we only consider finite horizon problem. It 
is described as follows. 

Let iV be a fixed integer and let c and (p be two measurable function, the optimal control 
problem in finite horizon is to consider what is called the "optimal cost" : 



minE 



'N-l 



k=0 



(61) 



If there is some strategy which realizes the minimum, this strategy is called the "optimal 
strategy". Let us investigate the classical result in stochastic control for the optimal 
strategy in this case. 

For this we define: 



min E 



.j=k I 



Remark The function c and are determined by the optimization constraint imposed 
by the experience. The equation which appears in the following theorem is called the cost 
equation and the function c and </> are called cost function. 

Theorem 9 Let U he a compact set and suppose that c is a continuous function. The 
solution of: 



V\p) 
V^{p) 

give the optimal cost: 

V\p) 



(62) 



minE 



'N-l 

^c{k,p],Uj) + (f){p'^, 

.j=k 



Pn = P 



The optimal strategy is given by: 



■■ P^ uUp) e argmin{p^+i(p)7^o' (p) + qk+i{p)T-C'^' (p) + c{k,p,Uk)}. (63) 
Furthermore this strategy is Markovian. 

Proof: The proof is based of what is called dynamic programming in stochastic control 
theory. Let u be any strategy and let V defined by the formula ( |62l) . we have 



mv'^' {Pk^i)-vHPk))/^{p:,^<k}] 
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then we have 



N-l 

Y.E[V'^\p-,_,,)-V\p-, 



k=0 
N-l 



N-l 

> — E [c{k, p", Uk)] (by definition of the min). 



k=0 



Hence for all strategy u, we have 



Vip) < E 



'N-l 



^c{k,p'^,Uk) + (j){pN) 



k=0 



Moreover we have equality if we choose the strategy given by the formula 
strategy is Markovian because the function c depends only on pk at time k. 



This 

□ 



The system (l62l) which describes the cost equation is called the discrete Hamilton- Jacobi 
Bellman equation. 

The fact that the optimal strategy is Markovian is another justification of the choice 
of such model of control for the discrete quantum trajectory. This theorem claims that we 
need just Markovian strategy in order to solve the "optimal control" problem. 

The next last section is devoted to the same investigation in continuous time models of 
quantum trajectories. 



3.2.2 The Continuous Case 

In the third section, we have proved the Poisson and the diffusion approximation in quan- 
tum measurement theory. We have the diffusive evolution equation 

Pt = po+ f L{s,p:,u{s,p:))ds+ fe{s,p:,u{s,p:))dw,, (64) 

Jo Jo 

and the jump-equation is 

Pt = Po+ i?(s-,p"_,M(s-,p"_))rfs (65) 

^0 

+ / / Qis, p'^^,u{s-, p'^_))lo<^<Tr[Jis-,u{s-,P^.)m.)]^idx,ds), (66) 

Jo Jr 

where the functions L, B, R and Q are defined in Section 2. 
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In this section, we consider the same problem of "optimal control" as in the discrete 
case. Let {fl, J-', J-'t, P) be a probability space where we consider the diffusive equation 



Pt 



and the jump-equation 



Po + 



t pt 

L{s,p'^,Us)ds+ / Q{s,p'^,Us)dWs 
Jo 



pt 



Po+ R{s-,p'^_,Us-)ds 
Jo 

+ 11 Q{s,p'^_,Us.)lo<x<Tr[J{s^,u,^){p^_)]N{dx,ds) 

Jo Jr 

where the strategy u = (ut) is just supposed to be a function J^t adapted (not only Marko- 
vian). In the case where J^t corresponds to the filtration generated by the process (pt), we 
recover the same definition as the discrete case. Concerning existence and uniqueness of a 
solution, with the condition (^31) of Section 2.1 for the functions L, R and 9 the previous 
equations admit a unique solution. Furthermore the solution takes values in the set of 
states on Ho- The set of all admissible strategy which satisfy the condition of adaptation 
is also denoted by U. The optimal control problem in this situation is expressed as follows. 

Let c and (p be two cost function. Let T > 0, the optimal problem in finite horizon is 
given by 



min E 



T 



cis,p'^,Us)ds + 0(p") 



As in the discrete model, we introduce the following function: 



V{t,p) = minE 



c(s, p", Us)ds + 0(p") / p^ = p 



(67) 



(68) 



which satisfies 



V{0,po) 



minE 



c{s,p'^,Us)ds + (t>{p'^) 



The function (!68l) represents the result of optimal control after t assuming pt = p. 

In this article, we just give the result for the optimal control problem for the diffusive 
case. A similar result for the Poisson case can be found in [TO] . 

As in the discrete case, it appears a continuous time version of the Hamilton- Jacobi- 
Bellmann Equation. The usual expression of this equation use the notion of infinitesimal 
generator of (p")- is described as follows in our context. A quantum trajectory (p") is 
considered as a process which takes values in with the identification of the state and 
the Bloch sphere Bi(R3) = {(x, y, z) e R^/x^ + y"^ + < 1}, that is. 



$ : Bi(M3) 
{x.y.z) 



M2(C) 
1 + X y + iz 
y — iz 1 — X 
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The map $ is injective and its range is the set of states. By considering that the functions L 
and G are apphcations from ]R_|_ x to M^, the stochastic differential equation concerning 
the diffusive case can be written as a system of stochastic differential equation on M'^ of 
the form: 

{Pt)i = Po+ f Li{s,p':,us)ds+ [ ei{s,p':,Us)dWs {1,2,3} 
Jo Jo 

where {pf)i (respectively 6j and Li) corresponds to the coordinate function of pf (respec- 
tively G and L). 

We introduce the 3x3 matrix 11 defined by Uij = QiQj. The infinitesimal generator 
of the process (p") acts on the functions / which are and bounded in the following 

way 

A-'^m = \ y: n^^(^'^'-)^ + (69) 

i,j=l * i=l * 

for all t > 0, M G M and x G M'^. In particular if m is a fixed constant, let (pt) be the 
solution of 

{pt)i = Po+ / Li{s,ps,u)ds + / <di{s,ps,u)dWs iG {1,2,3}. 
Jo Jo 

Hence for all function / which is and bounded, the following process 

Mt = f{pt)-f{po)- f A^^'f{ps)ds 

Jo 

is a martingale for the filtration generated by (pt). 

The following theorem express the result in optimal control for the diffusive quantum 
trajectory. 

Theorem 10 Suppose there is a function {t,p) V{t,p) which is in t and in p 
such that: 

^^ + min„e«{^"-*F(t,p)+c(t,p,M)} = . . 

V{T,p) = <j>{p) ^^^^ 

where A^'^fi^x) is defined by the expression (fgPj) . The function V gives the solution of the 
optimal problem, that is, 



V{t,p) = minE 



c(s, p^, Us)ds + (f){p^) / = p 



Furthermore if the strategy u defined by 

u\t,p) G argmin{^"'V(t,p) +c(t,p,M)} (71) 

is an admissible strategy then it defines an optimal strategy. Moreover this strategy is 
Markovian. 
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The equation fl70p is the Hamilton- Jacobi-Bellmann equation in the continuous case. 

A proof of this theorem can be found in [22] or [2H]- The interest of such theorem in 
our context is to notice that the optimal strategy is Markovian, this confirms the choice of 
such strategy in the model of quantum trajectories with control. 

A similar result holds for the Poisson case. The infinitesimal generator for such process 
is given in [13], explicit computations can be found in [27]. 



References 

S. Attal. Approximating the Fock space with the toy Fock space. In Seminaire de 
Probabilites, XXXVI, volume 1801 of Lecture Notes in Math., pages 477-491. Springer, 
Berlin, 2003. 

S. Attal and A. Joye. The Langevin equation for a quantum heat bath. J. Funct. 
Anal, 247(2):253-288, 2007. 

S. Attal, A. Joye, and C.-A. Fillet, editors. Open quantum systems. Ill, volume 1882 
of Lecture Notes in Mathematics. Springer- Verlag, Berlin, 2006. Recent developments. 
Lecture notes from the Summer School held in Grenoble, June 16-July 4, 2003. 

S. Attal and Y. Fautrat. From repeated to continuous quantum interactions. Ann. 
Henri Poincare, 7(1):59-104, 2006. 

S. Attal and C. Fellegrini. Contributions to quantum trajectory, preprint, 2007. 

A. Barchielli and G. Lupieri. Quantum stochastic models of two-level atoms and 
electromagnetic cross sections. J. Math. Phys., 41(ll):7181-7205, 2000. 

A. Barchielli and G. Lupieri. Instrumental processes, entropies, information in quan- 
tum continual measurements. Quantum Inf. Comput., 4(6-7) :437-449, 2004. 

A. Barchielli, A. M. Faganoni, and F. Zucca. On stochastic differential equations and 
semigroups of probability operators in quantum probability. Stochastic Process. AppL, 
73(l):69-86, 1998. 

L. Bouten, M. Guta, and H. Maassen. Stochastic Schrodinger equations. J. Phys. A, 
37(9):3189-3209, 2004. 

L. Bouten and R. van Handel. Quantum filtering: a reference probability approach, 
2005. 

L. Bouten, R. van Handel, and M. James. An introduction to quantum filtering, 2006. 

N. Bruti-Liberati and E. Flaten. On the strong approximation of jump- 
diffusion processes. Research Faper Series 157, Quantitative Finance Re- 
search Centre, University of Technology, Sydney, April 2005. available at 



http :/ /ideas . repec . org /p/uts/rpaper / 1 5 7 . html 



36 



[13] P. Cheridito, D. Filipovic, and M. Yor. Equivalent and absolutely continuous measure 
changes for jump-diffusion processes. Annals of Applied ProhaBILITY, 15:1713, 2005. 

[14] E. B. Davies. Quantum theory of open systems. Academic Press [Harcourt Brace 
Jovanovich Publishers], London, 1976. 

[15] S. Haroche and J.-M. Raimond. Exploring the quantum. Oxford Graduate Texts. 
Oxford University Press, Oxford, 2006. Atoms, cavities and photons. 

[16] J. Jacod and P. Protter. Quelques remarques sur un nouveau type d'equations 
differentielles stochastiqucs. In Seminar on Probability, XVI, volume 920 of Lecture 
Notes in Math., pages 447-458. Springer, Berlin, 1982. 

[17] J. Jacod and A. N. Shiryacv. Limit theorems for stochastic processes, volume 288 of 
Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathe- 
matical Sciences]. Springer- Verlag, Berlin, second edition, 2003. 

[18] R. Kadison and J. R. Ringrose. Fundamentals of the theory of operator algebras. Vol. 

I, volume 15 of Graduate Studies in Mathematics. American Mathematical Society, 
Providence, Rl, 1997. Elementary theory. Reprint of the 1983 original. 

[19] R. Kadison and J. R. Ringrose. Fundamentals of the theory of operator algebras. Vol. 

II, volume 16 of Graduate Studies in Mathematics. American Mathematical Society, 
Providence, RI, 1997. Advanced theory. Corrected reprint of the 1986 original. 

[20] T. G. Kurtz and P. Protter. Weak limit theorems for stochastic integrals and stochastic 
differential equations. Ann. Probab., 19(3):1035-1070, 1991. 

[21] T. G. Kurtz and P. Protter. Wong-Zakai corrections, random evolutions, and sim- 
ulation schemes for SDEs. In Stochastic analysis, pages 331-346. Academic Press, 
Boston, MA, 1991. 

[22] H. Kushner. Introduction to stochastic control. Holt, Rinehart and Winston, Inc., 
New York, 1971. 

[23] K. R. Parthasarathy. An introduction to quantum stochastic calculus, volume 85 of 
Monographs in Mathematics. Birkhauser Verlag, Basel, 1992. 

[24] Y. Pautrat. Prom Pauli matrices to quantum Ito formula. Math. Phys. Anal. Geom., 
8(2):121-155, 2005. 

[25] C. Pellegrini. Existence, uniqueness and approximation of stochastic Schrddinger equa- 
tion: the Poisson case, submitted, 2007. 

[26] C. Pellegrini. Existence, Uniqueness and Approximation of a Stochastic Schrddinger 
Equation: the Diffusive Case, to appear in The Annals of Probability, 2008. 



37 



[27] C. Pellegrini. Markov Chains Approximation of Jump-Diffusion Quantum Trajectories. 
preprint, 2008. 

[28] H. Pham. On some recent aspects of stochastic control and their applications. Prob- 
ability Surveys, 2:506, 2005. 

[29] P. E. Protter. Stochastic integration and differential equations, volume 21 of Appli- 
cations of Mathematics (New York). Springer- Verlag, Berlin, second edition, 2004. 
Stochastic ModeUing and Applied Probability. 

[30] R. van Handel, J. K. Stockton, and H. Mabuchi. Feedback control of quantum state 
reduction. IEEE Transactions on Automatic Control, 50:768, 2005. 

[31] Didier Dacunha-Castelle and Marie Dufio. Probability and statistics. Vol. II. Springer- 
Verlag, New York, 1986. Translated from the French by David McHale. 



38 



